Quasiparticle Density of States of Clean and Dirty s-Wave Superconductors 

in the Vortex State 
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The quasiparticle density of states (DOS) in the vortex state has been probed by specific heat mea- 
surements under magnetic fields {H) for clean and dirty s-wave superconductors, Y(Nii_a;Pti:)2B2C 
and Nbi_a;Tai:Se2. We find that the quasiparticle DOS per vortex is appreciably //-dependent 
in the clean-limit superconductors, while it is //-independent in the dirty superconductors as ex- 
pected from a conventional rigid normal electron core picture. We discuss possible origins for our 
observations in terms of the shrinking of the vortex core radius with increasing H. 



PACS numbers; 74.25.Bt, 74.60.Ec, 74.70.Ad 

Determination of the electronic structure of magnetic 
vortex lines in type-II superconductors has been an issue 
of long standing interest. Recent attention has focussed 
on its relation with the symmetry of the Cooper pairs: 
For s-wave superconductors with an isotropic gap, there 
exist low-lying bound quasiparticles around the center of 
the vortex lines, where the superconducting gap vanishes, 
and the vortex line has been simply viewed as a "rigid" 
normal electron cylinder with a radius of the coherence 
length £_. ||] On the other hand, for d-wave superconduc- 
tors, quasiparticles are delocalized and spread outside 
the vortex cores because of the presence of the gapless 
region on the Fermi surface. This difference manifests 
itself in, for example, a magnetic field (i/) dependence 
of the quasiparticle density of states (DOS) at the Fermi 
level No{H); No{H) ex H for s-wave superconductors, 
while A^o (H) ex vi? for d-wave superconductors with line 
nodes. [0 However, even in s-wave superconductors, the 
following unsolved issues may potentially alter the above 
simple picture of the vortex lines: Shrinking of the core 
radius with decreasing temperature (T) , due to the ther- 
mal population of the quantized levels for the bound ex- 
citations, has been predicted theoretically by Pesch and 
Kramer. ||^ Shrinking of the core radius with increasing 
H has been proposed based on scanning tunneling spec- 
troscopy (STS) y and muon spin rotation (^SR) Q mea- 
surements, though the physics behind this has not been 
fully explored yet. The de Haas - van Alphen (dHvA) 
effect has been observed in the mixed state where the 
cyclotron radius is much larger than ^. ||,0] This sug- 
gests the existence of delocalized quasiparticles even in 
the s-wave superconductors under high magnetic fields. 

In this letter, we report the low-temperature spe- 
cific heat of pure and alloyed Y(Nii_a;Pta;)2B2C and 
Nbi_a:Taa;Se2 single crystals under magnetic fields, fo- 
cusing on the H dependence of the quasiparticle DOS. 
We will show that the H dependence of the quasiparti- 
cle DOS at the Fermi level No{H) is significantly influ- 
enced by the presence of disorder: While the DOS shows 
the No{H) ex H behavior in the alloyed samples, No{H) 



shows significant deviation from the Nq{H) cx H behav- 
ior in the pure samples. We will discuss possible mecha- 
nisms of the deviations from the //-linear behavior and 
their relation with quasiparticle scattering. 

Single crystals of Y(Nii_a;Pta;)2B2C were grown 
by a floating zone method, [gj while iodine- assisted 
chemical transport was employed for the growth of 
Nbi_:rTaa;Se2. |y] Table I lists the superconducting and 
normal state parameters for the samples used in this 
study. As evident from the table, the pure samples, 
YNi2B2C and NbSe2, are in the clean limit (^ ^ /, 
where I is the mean free path), while the alloyed sam- 
ples, Y(Nio.8Pto. 2)2620 and Nbo.8Tao.2Se2, are in the 
intermediate range near the dirty limit (C ^ 0- To deter- 
mine Nq(H) in magnetic fields, we have measured low- 
temperature specific heat (C) from 1.4 K to 20 K under 
magnetic fields up to 12 T using a thermal relaxation 
calorimeter, mn The relative resolution of the measure- 
ments was better than 0.5 %, and the absolute accuracy 
determined from the measurement of a Cu standard was 
better than 1 %. 

Specific heat data for Y(Nii_a;Pta;)2B2C under various 
magnetic fields parallel to the c axis are shown in Fig. 
1, plotted as C/T versus T^ for two samples (a) x — 
0.0 and (b) x = 0.2. The measurements were performed 
with increasing temperature after field cooling from a 
temperature well above the superconducting transition 
temperature Tc- For the pure sample, a clear jump is 
seen in the zero-field data at Tc = 15.4 K. The upper 
critical field Hc2{T) was determined from the specific 
heat data for a fixed H. The value at T = 0, /Jc2{0), 
was estimated by fitting the low-temperature data be- 
low 5 K to Hc2{T) = /?c2(0)[l - h{T/Tc)% where h is a 
constant, giving an estimate of i/c2(0) — 8.0 T. This is 
consistent with a reported value of i/c2(0) = 8.1 T. |l|] 
The least squares fit of the 9 T data, which represent the 
normal- state specific heat, to C = jt^T -\- [3T^ between 
1.4 K and 5 K provides a Sommerfeld constant 7n = 20.6 
mJ/K^mol. These values agree reasonably with those re- 
ported previously. |13] For the x — 0.2 sample, we obtain 



Tc = 12.1 K, 7n = 14.7 mJ/K^mol, and Fc2(0) = 4.3 T. 

Under magnetic fields, a finite amount of T-linear 
term, j{H)T, appears in the specific heat well below Tc as 
shown in Fig. 1. The 7(i?) values at each field were esti- 
mated by extrapolating the C/T versus T^ data between 
1.4 and 2.5 K linearly to T = 0. The zero-field value 
7(0) for the pure sample is negligibly small, indicative of 
the high quality of the sample. A small residual 7(0) of 
about 0.4 mJ/K^mol is seen for the x — 0.2 sample, which 
very likely originates from a tiny amount of normal-state 
secondary phase. In order to see the magnetic field de- 
pendence, 7(i/)/7N is plotted in Fig. 3(a) as a function 
of H/Hc2- A substantial difference between the two sam- 
ples is evident. For the pure sample, j{H) first shows a 
marked increase at low fields, then gradually approaches 
the normal state value. Notably, j{H) reaches a value 
almost half that of 7n at 1 T (7J/7?c2 ^ 0.13). We ob- 
served essentially the same 7?-dependence oi j{H) in the 
field direction H J- c. Analogous nonlinear dependence 
of 7(iJ) was also observed for LuNi2B2C. |1^ In contrast, 
for the X = 0.2 sample, 7(il) increases linearly with H 
as 'j{H) = 7(0) -t-7Niif/-ffc2 over a wide field range below 
Hc2- l{H) does not show any field dependence above 
Hc2, indicating that this region is indeed in the normal 
state. 

Qualitatively similar behavior for "f{H) was observed 
in Nbi_a;Taa;Se2. Shown in Fig. 2 are the specific heat 
data for Nbi_a;Ta2;Se2 under magnetic fields applied par- 
allel to the c axis. For the pure sample, x — 0.0, we 
obtained T^ = 7.3 K, 7n = 19.3 mJ/K^mol, and -ffc2(0) 
= 4.4 T. For the alloyed sample, x — 0.2, we obtained 
Tc = 5.1 K, 7N = 15.1 mJ/K^mol, and i?c2(0) = 3.2 T. 
The magnetic field dependence of 7(i?) was estimated in 
the same manner as mentioned above, and is displayed in 
Fig. 3(b). A nonlinear behavior oi j{H) is noticeable for 
the pure NbSe2 sample, though not as significant as that 
observed in YNi2B2C. Similar nonlinear behavior oi^ {H) 
in the pure sample was reported by Sanchez et al. [ |l3yi4| 
In contrast, for the alloyed sample, Nbo.8Tao.2Se2, ^{H) 
exhibits a linear dependence on H below iJc2 as observed 
in Y(Nio.8Pto.2)2B2C. 

Since the number of vortices is scaled by H, the nonlin- 
ear behavior of 7(7?) in the pure superconductors implies 
that the DOS per single vortex depends on H when the 
sample is clean. This contradicts the conventional pic- 
ture that the vortex line should be viewed as a "rigid" 
normal electron cylinder, where the bound quasiparti- 
cles give rise to a if -independent DOS proportional to 
iVoTT^ ^ per vortex line, where A^o is the DOS in the nor- 
mal state. However, once the electron mean free path I 
is substantially reduced due to strong impurity scatter- 
ing and / < ^ < intervortex distance i?, the conventional 
picture appears to hold true. We point out that the un- 
usual 7(if ) behavior and the strong infiuence of disorder 
are universally observed for s-wave superconductors, al- 
though the detailed 7(il) data are available only for a 



limited number of compounds: A nonlinear dependence 
of 7(i?) was reported on a clean s-wave superconductor 
CeRu2, ||l5[ where the deviation is compatible with that 
observed in NbSe2. By contrast, a linear behavior was re- 
ported on a dirty superconductor Nb77Zr23. 10 Both H- 
linear and nonlinear ^{H) were observed on VaSi, [^ 17 
which is very likely due to different degrees of disorder in 
the samples. 

We suggest that the apparent breakdown of the "rigid" 
normal electron picture originates from the H depen- 
dent core radius. In the conventional picture, the vor- 
tex core is assumed to be rigid. However, if the core 
radius is H dependent, the quasiparticle DOS per vor- 
tex should also be H dependent, resulting in a nonlinear 
j{H). Indeed, direct measurement of the local DOS in 
NbSe2 by STS gives evidence for the shrinking of the 
core radius. Q Using the experimentally observed 7n 
and 7(-ff), we estimate the effective core radius py by 
defining 2TTpy{H) 7n — j{H)/{H/^o), where $0 is the 
flux quantum. Note that pv(^c2) = C(0)- The variation 
of pv (H) estimated in this way is shown in the inset to 
Fig. 3, and compared with those determined by STS for 
NbSe2. The H dependence of py from these two in- 
dependent probes appears to agree reasonably with each 
other. The difference in the magnitudes should not be 
taken as significant since there is an arbitrariness of a 
factor in the definition of py. 

The shrinking of the core radius may be supported by 
an anomalous temperature dependence of Hc2iT) near 
Tc. In Fig. 4(a), Hc2{T) curves deduced from the spe- 
cific heat data are shown for the pure YNi2B2C where 
the deviation of j{H) from the linear behavior is most 
significant. A pronounced positive curvature is clearly 
seen near Tc. This may be understood in terms of the 
shrinking of py (and hence ^), which would lead to 
an enhancement of -ffc2(T) with increasing H through 
the relation Hc2 = ^o/27r^^. Indeed, from the sup- 
pressed initial slope at Tc, diJc2/dT|y^ = -0.30 T/K, 
we can obtain an estimate for ^(0) ~ 115 A using 
C(0) = 0.54[-$o/rc(diic2/dT|j,J]i/^ which is substan- 
tially larger than ^(0) ~ 65 A, but close to the estimate 
of the enhanced py ~ 130 A from j{H) at low mag- 
netic fields limit (Fig. 3). Furthermore, as clearly seen in 
Fig. 4(a), the alloyed sample Y(Nii_a;Pta;)2B2C does not 
show any noticeable positive curvature, consistent with 
the iJ-independent py . |1^ A small positive curvature of 
Hc2{T) near Tc can be seen in the pure NbSe2 as shown 
in Fig. 5(b) (See also ref. (l|) and CeRu2, ^ which 
is consistent with the observation of a definite but less 
pronounced deviation from the H-linear "f{H) in these 
materials. 

By introducing a certain (repulsive) vortex-vortex in- 
teraction in the vortex lattice states, it was shown theo- 
retically that the core radius indeed shrinks with increas- 
ing H. I4l@,n3] The physical origin of the vortex interac- 



tion is not clear yet. We speculate that this interaction 
is mediated by a coherent transfer of quasiparticles be- 
tween the cores. The observation of dHvA oscillations 
in the mixed state of clean YNi2B2C @] and NbSe2, ||] 
where the cyclotron radius is by far larger than the core 
radius, clearly indicates the presence of delocalized and 
highly coherent quasiparticles outside the vortex cores. 
These delocalized quasiparticles may give rise to a cer- 
tain vortex-vortex interaction and result in a shrinking 
of the core radius with increasing H . Then the critical 
sensitivity of the ^{H) behavior can be naturally under- 
stood, since the coherent motion of these quasiparticles 
will be substantially suppressed by the impurity scatter- 
ings, when I ^ S, < R- One may argue that the DOS asso- 
ciated with the delocalized quasiparticle itself may be the 
origin of nonlinear ^{H) behavior. Hedo et al. recently 
proposed that the nonlinear behavior of ^(H) in CeRu2 
originates from the contribution of delocalized quasipar- 
ticles with momentum k perpendicular to H. [|l5| The 
discussion is based on the theoretical study by Brandt 
et al., Ig^ who showed that near iJc2, while excitations 
with k parallel to H have a gap, those in planes per- 
pendicular to H are gapless due to the formation of the 
vortex line lattice. However, the theory by Brandt can 
be applicable only near -ffc2, and another theoretical ap- 
proach with emphasis on the low field regime indicates 
that the local DOS outside the core is almost zero when 
H/H,2 < 0.4. ^ 

Finally, we would like to comment on the pronounced 
nonlinear 7(i/) behavior in pure YNi2B2C compared 
with that of NbSe2. We suspect that the presence of 
substantial gap anisotropy plays a role in YNi2B2C. The 
presence of a reduced gap region is suggested by the 
dHvA measurement. H As seen in Fig. 1(a), the zero- 
field specific heat of pure YNi2B2C shows T^ behavior 
rather than a thermally activated behavior down to low 
temperatures. This also supports the presence of a sub- 
stantially small gap region on the Fermi surface, although 
d- or p- wave pairing is unlikely in YNi2B2C because of the 
insensitivity of the superconductivity to nonmagnetic Pt 
substitution. The reduced gap region on the Fermi sur- 
face may enhance quasiparticle transfer between the vor- 
tices and may result in the pronounced nonlinear ^{H) in 
terms of the above speculative picture. In NbSe2, zero- 
field specific heat exhibits thermally activated behavior 
at low temperatures, indicative the presence of a rela- 
tively isotropic gap. This is consistent with the smaller 
deviation from the linear behavior of ^{H) and positive 
curvature oi Hc2{T) in this compound. 

In summary, we have measured the specific heat of 
Y(Nii_2;Pta;)2B2C and Nbi_2:Taa;Se2 in the vortex state. 
Our results show that, in clean superconductors, the co- 
efficient of the T-linear term in the specific heat, 7(-ff), 
shows nonlinear behavior with H, suggesting that the 
DOS per vortex is H dependent, while in dirty supercon- 
ductors, j{H) is linear, suggesting a constant DOS per 



vortex. We suggest that these results may be understood 
in terms of the shrinking of the vortex core radius with 
increasing magnetic field. To strengthen this scenario, it 
is critically important to examine the H dependence of 
pv in dirty superconductors by direct probes such as STS 
and /iSR. Such attempts are now in progress. 
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FIG. 1. Specific heat divided by temperature C /T for 
(a) YNiaBaC and (b) Y(Nio.8Pto.2)2B2C, as a function of T^ 
under various magnetic fields. The solid lines are guides to 
the eye. 
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FIG. 2. Specific heat divided by temperature C /T for 
(a) NbSe2 and (b) Nbo.8Tao.2Se2, as a function of T^ under 
various magnetic fields. The solid lines are guides to the eye. 
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FIG. 3. The coefficient of T-finear specific heat at 

low temperatures normalized by the Sommerfeld constant, 
7(if)/7N, as a function of reduced magnetic fields, H/Hc2, 
for (a) Y(Nii_^Pt^)2B2C and (b) Nbi_:,Ta^Se2. The solid 
lines are guides to the eye. The insets show the ef- 
fective vortex core radius, pv, determined by assuming 
27rpv(^)^7N = 'y{H)/{H/^o), as a function of H/Hc2- pv 
of NbSe2 determined by STS || at T/Tc = 0.33 are shown 
(open squares) for comparison. 
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FIG. 4. Temperature dependence of the upper critical 
field i/c2(T) for Y(Nii_^Pt^)2B2C and Nbi_^Ta^Se2. The 
solid lines are guides to the eye. The broken lines indicates 
the fit of the initial slope of Hc2{T) at Tc- 



I. TABLE 



TABLE L Superconducting and normal state parameters 
for the Y(Nii_i:Pta;)2B2C and Nbi_a;Taa;Se2 samples studied. 
The mean free path I was calculated from the resistivity po 
just above Tc by using / = h{S'K^) /e^n^'^po- We used an 
electron density of 3 x 10^^ e/cm^ for Y(Nii_j;Pt2;)2B2C and 
1.6 X 10^^ e/cm^ for Nbi-^Taa;Se2. The Ginzburg Landau 
coherence length parallel to the ab plane, ^{0), was determined 
fromHc2(0). 
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